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Abstract 

Training and recovery of exchange bias in FeNi/Cu/Co/FeMn spin valves 
have been studied by magnetoresistance curves with field sweep rates from 
1000 to 4800 Oe/s. It is found that training and recovery of exchange field 
are proportional to the logarithm of the training cycles and recovery time, 
respectively. These behaviours are explained within the model based on 
thermal activation. For the field sweep rates of 1000, 2000 and 4000 Oe/s, 
the relaxation time of antiferromagnet spins are 61.4, 27.6, and 11.5 in the 
unit of ms respectively, much shorter than the long relaxation time (~ 10 2 s) 
in conventional magnetometry measurements. 
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1. Introduction 

The exchange bias (EB) effect in ferromagnetic/ antiferromagnetic sys- 
tems have been intensely studied in the last decade because of their physical 
complexity and important applications [H, 0]. The technological importance 
lies in the pinning effect of the antiferromagnet (AFM) layers in which the 
hysteresis loop of the ferromagnet (FM) can be shifted away from the origin 
point by the amount of the exchange field (He), and is usually accompanied 
with an enhanced coercivity (He)- Changes of and Hq are accordingly 
directly related to the spin configuration of the AFM layer through the ex- 
change coupling [3J. Among the variety of effects related to the EB phe- 
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nomenon, the training effect is an important effect that reflects the AFM 
spin dynamic process during repeated hysteresis loops. It is ascribed to that 
the spin structure of the AFM layer deviates from its equilibrium configura- 
tion and approaches another equilibrium triggered by subsequent reversals of 
the FM magnetization. Nowadays, studies of AFM spin dynamic behaviours 
with trai ning effect in both experiments and theories have been widely re- 
ported p-13]. Because most of studies are limited to long timescales (> Is), 
by the usually quite long measurement time in magnetometry approaches, 
the relaxation time of AFM spin are usually reported in second timescale 
(~ 10 2 — 10 4 s)[3-@]. In contrast, at shorter measurement timescales the re- 
laxation time of exchan ge b ias system was demonstrated to cover a wide 



range ( ~ 10 8 — 10 11 s) [15l-ll8|. which has been ascribed to the magnetiza- 



tion reversal mechanism of FM layer |14j . Hence, the report only on AFM 



spin dynamic behaviour in the millisecond timescale is still sparse. In ad- 
dition, recently attempt frequencies up to 10 12 Hz in AFM layer have been 



reported [19], which indicated the much shorter relaxation timescale of AFM 
spin than earlier anticipated. Therefore, it is necessary and interesting to 
study the AFM spin dynamic process at short timescale (technologic impor- 
tance < Is). 

In this work, we have studied the EB training and recovery behaviours 
at the millisecond timescale based on the electrical transport measurements 
in FeNi/Cu/Co/FeMn spin valves. The experiments show that at high field 
sweep rates recovery time of exchange field after training procedures is three 
orders of magnitude shorter than the values observed by usual magnetometry 
techniques, and the relaxation of magnetoresistance (MR) is demonstrated 
in the millisecond timescale. These clearly indicate that AFM spin dynamic 
behaviours can be studied and resolved down to the millisecond timescale 
utilizing the ordinary resistance measurements. 



2. Experiment and Results 

The spin valves of Si (001)/Cu (10 nm)/Fe 20 Ni 80 (3 nm)/Cu (3 nm)/Co 
(3 nm)/FeMn (8 nm)/Ta (3 nm) were prepared by a magnetron sputtering 
system. The base pressure was 2xl0 -5 Pa and the Ar pressure was 0.3 Pa 
during the deposition. The 10 nm Cu buffer layer was used to stimulate 
the fee (111) preferred growth of the FeMn layer in order to enhance the 
EB. A magnetic field of 130 Oe was applied in the film plane during depo- 



2 



sition to induce the uniaxial anisotropy and thus the EB. Magnetoresistance 
(MR) measurements were performed to probe the switching behaviours of 
the pinned layer for different subsequent hysteresis loops. The magnetic field 
was provided by home-built Helmholtz coils, and MR was measured in real- 
time system with 2 M/s sampling rate. To study training and recovery of 
the EB, we first performed forty consecutive MR measurements with a fixed 
field sweep rate to characterize the training procedures. Then we stopped 
the magnetic field sweep with an waiting time t. Finally ten consecutive 
MR measurements with the same field sweep rate were measured in order to 
observe and confirm the EB recovery. For each sweep rate, t varied from 0.1 
to 10 s. 

The spin valves MR curves of the training and recovery effects at the 
1 st , 40 th and 41 s * cycles with the field sweep rate of 4000 Oe/s are displayed 
in Fig. [TJ (a). At large negative field the Co and FeNi magnetizations are 
parallel and pointing down. When the field is increased above the switch 
field of the Co layer, about -110 Oe, the Co magnetization reverses and re- 
sistance switches from low value (-1) to high value (+1). When the field is 
further increased above the switch field of the FeNi layer, about -15 Oe, its 
magnetization reverses, the two magnetizations become parallel once more 
but this time pointing up, and resistance switches to low value (-1). If the 
field is then decreased, the two magnetizations will remain parallel until the 
negative switch field of the FeNi layer is reached at -25 Oe, when its mag- 
netization reverses and resistance switches to high value. When the field is 
further reduced and reverses the Co magnetization, the two magnetizations 
align in parallel, and resistance changes to its low value. For all MR curves 
the hysteresis loops of the Co layer are shifted and fully separated from the 
hysteresis loops of FeNi layer due to the FeMn pinning effect, therefore the 
MR curves directly reflect the switching behaviours of the Co and the FeNi 
layers in detail [20]. Comparing the hysteresis loops of the Co layer in the 
I s ', and 40 th MR curves, the switching field of the descent branch shifts more 
sharply than that of the ascent one, demonstrating the asymmetric magne- 
tization reversal. However, after the magnetic field sweep is stopped for 1 s, 
a recovery is observed in the 41 st MR curve. It contrasts to the behaviour 
in the case of normally low field sweep rate, in which substantial recovery 
was only observed after several hours of waiting time (3). The is plotted 
as function of cycles n in Fig. QJb). The gradually decreases with the 
cycle n, has an obvious resilience after 1 s waiting time and finally decreases. 
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For the training procedure, the versus n is fitted by a linear functions of 
l/y/n , e~ a06ri and ln(n). It is found that the logarithm function yields the 
best fit, except for initial point n—1 (8l. [Toj] . 

To further study the recovery of the trained EB, we measured the re- 
covery rate R as a function of t at different field sweep rates, where R = 
[F E (41) - tf E (40)]/[F E (l) - # E (40)] x 100(%). Figure Ha) shows the de- 
pendence of R on t at different field sweep rates. The R increases with 
the increasing t as a linear function of log(t). More remarkably, for a fixed 
waiting time t, R correspondingly increases with the increasing field sweep 
rate. This logarithm behaviour is in a good agreement with the previous 
experiments in NiFe/FeMn system, while the recovery rate is several orders 
of magnitude faster than the value in the low field sweep rate [4]. The slope 
and the intercept as a function of the field sweep rate are shown in Fig. |2](b). 
The slope displays little change with different field sweep rates whereas the 
intercept increases greatly as the field sweep rate increases in approximate 
linear function of the logarithm of the field sweep rate. 

To investigate the dynamic behaviour of the EB with high resolution, 
we observed the evolution of MR after setting the magnetic field from the 
positive saturation field to -210 Oe (the point A in Fig. [U^a)) near the switch 
field. As shown in Fig. [3j MR initially decreases sharply and then gradu- 
ally reaches a constant. The small fluctuations in the curves are caused by 
50 Hz AC noise in the amplifying circuit. Remarkably, a crossover of the 
normalized MR from positive to negative has been observed, demonstrating 
the reversal of the magnetization of the Co layer. It is possible to link the 
time dependence of MR with the magnetic viscosity in the Co/FeMn bilayers 
[t| , in which the magnetization of the pinned layer gradually reverses due to 
the thermally activated process in Co/FeMn bilayers. Because the reversal 
process in the EB at the first cycle consists in the single domain wall motion 
[6], the change of MR here is proportional to the amounts of the reversal 
magnetization in the pinned layer. Shown as the solid line in Fig. [3j the 
evolutions of MR are described well by a first order exponential decay. From 
fitting the data, we extracted the relaxation times r are 11.5, 27.6, and 61.4 
ms for the field sweep rate at 4000, 2000 and 1000 Oe/s respectively. This 
is again in contrast to the long relaxation time (~800 s) in the conventional 
approaches H. One can also note the relaxation time decreases with the 
increasing field sweep rate. 
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3. Discussion 



The above results show that the recovery and relaxation of the EB at 
high field sweep rates are faster than that earlier observed 0-0|- Below we 
will interpret the experimental results in conventional models for AFM and 
training effects. 

Firstly we consider the change and magnitude of the relaxation time 
constants at different field sweep rates shown in Fig. [3j The time con- 
stant for the relaxation can be described by an ordinary Arrhenius law 
t = v^ 1 exp(£^ (T //c B T), where v a is the attempting frequency and E a = KV 
represents the AFM energy barrier, K is the AFM anisotropy and V is the 
AFM grain volume. According to the AFM grain volumes distribution, we 
can divide the E a into three different categories [lOj]: i) small E a (small grain 
size), which follows the FM magnetization at the timescale of the experiment. 
ii) medium energy E a (medium grain size) which will determine the EB dy- 
namics at the timescale we investigate. Hi) large E a (large grain size), which 
is a stable configuration over the timescale of the experiment. Assuming a 
typical uniaxial anisotropy constant of 1 x 10 6 erg/cm 3 and v a to be 1 x 10 9 
Hz, then the average grain size of category (ii) is correspondingly about 9 
nm extension, based on the relaxation time in Fig. [3j Accordingly, that re- 
laxation time decreasing with the increasing field sweep rate, demonstrates 
an apparent increase in attempt frequency v a . 

The EB recovery and relaxation at high field sweep rates can still be ex- 
plained well with the model based on thermal activation 0, [23[ . As shown 
in Fig. 12(a), the logarithm time recovery relationship indicates a thermally 
activated reversal process involving the AFM spin configuration. To explain 
our data, the activation energy spectrum model simply based on a two-level 



system is adopted [23|, |24[. In our case the two level system represents an 
individual AFM grain or domain switching from a positive to a negative ex- 
change energy with respect to the FM layer. For the system with a wide 
energy barrier distribution, AH^ can be expressed in terms of the AFM acti- 
vation energy spectrum q(E): AH^ = q(E)kBTln(v (T t) , which is taken from 
equation (1) in ref Q. According to the equation, the slope observed for all 
field rates in Fig. [2(b) is a constant due to the same q(E), while the intercept 
variation is mainly due to the different activated AFM energy ranges and the 
time delay at the different field sweep rates. 

Finally, for the training process is proportional to ln(n) at high field 
sweep rates in Fig. 1(b), which can be compared to the usual power law 



5 



1/y/n) and the exponential (e an ) relationships. We model this through 



following Binek et al [21|, |22j. At beginning, the equilibrium AFM interface 
magnetization is defined S AFM = lim^oo Safm (n). Each positive and neg- 
ative deviation SS n = Safat fa) — S AFM of the AFM interface magnetization 
from its equilibrium value will increase the total free energy F of the sys- 
tem by AF. The relaxation of the system towards equilibrium is determined 
by the Landau-Khalatnikov (LK) equation 25]: ^Safm = —OAF /8Safmi 
where £ is a phenomenological damping constant and AF is the function of 
5S. In Binek's model under the assumption AF(SS) = AF(—5S) , a se- 
ries expansion of AF up to the fourth order in 5S yields AF = ^a(5S) 2 + 
^b(5S) 4 + 0(5S) 6 . Evaluating the free energy expression with a leading term 



of 2 nd and 4 th order in SS will result in the e~ an 22| and l/\/n 21] evolution, 
respectively. 

However, our understanding of the system is that we have a nonvanishing 
odd order term. This is an effect of working at a time scale where we also have 
substantial coupling at the FM/AFM interface due to large grains that are 
too large to follow the oscillating exchange coupling of the FM. Instead that 
portion of the ensemble of grain will orient itself gradually according to the 
mean coupling induced by the FM in a monotonic fashion. Accordingly we 
also have to considered the expansion of AF also from first order of SS. We 
then assume SS: AF = /(ra)((5S') 1 + 0(SS) 2 , where the /fa) indicates that 
the change in the AFM interface magnetization and SS n is dependent on the 
training procedures n. By replacing S with [S(n+ 1) — S(n)]/r, with r being 
the relevant experimental time constant and the free energy expression of the 
first order into the LK equation, we obtain an implicit sequence equation: 
£'(S(n + 1) — S(n)) = —fin), where £' = £/r. The sum of this equation 
yields H E (n + 1) oc S(n + 1) = - £ /(n)/£'. 

Since we do not know the exact energy distribution of our system, we 
make a first order approximation assuming a constant distribution of areas, 
which leads to a linearly increasing activation energy with time, then the 
approximate development f(n) will follow a 1/n dependance. Using this ap- 
proximation, we obtain a logarithmic relaxation of the system that fits well 
with the observed dependence. Note, that already at n=6 the approxima- 
tion holds to within an error of less than 4%. We note that the training 
process H E is proportional to ln{n) has also been reported at low field sweep 
rates [8| , where the training speed is several orders of magnitude slower than 
the values reported here. 

In summary, for the AFM spins the relaxation time in the millisecond 
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timescale is demonstrated when the bilayers are exposed to high field sweep 
rates. This behaviour can be well explained in terms of a time constrained 
thermal activation. Our finding gives a new insight into the dynamic be- 
haviour of the AFM spins. 
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Figure 1: (a) The magnetoresistance curves used to map the training effect for 
FeNi/Cu/Co/FeMn spin vaive at the I s *, 40 th and 41 st (after 1 s waiting time) cycles 
with the held sweep rate of 4000 Oe/s. T@e resistance is dependent on corresponding 
magnetization configurations of FeNi (black arrow) and Co (red arrow) . (b) The exchange 
field He as a function of the number of cycles n. The blue dot, green dash dot and black 
dash lines are the fitted data with the 1/y/n , e _0 06n and ln(n), respectively. 




Figure 2: (a) The recovery of He as a function of the waiting time t with different sweep 
rates. The solid lines display the linear fits of the ln(t). (b) The slope and the offset values 
as a function of the field sweep rate. The solid lines are the linear fits of the logarithm of 
the field sweep rate. 
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Figure 3: The time dependence of the resistance after the external magnetic field is swept 
to -210 Oe (point A in Fig. Qla)) from positive saturation field with different field sweep 
rates (a) 4000 Oe/s (b) 2000 Oe/s (c) 1000 Qe/s. The solid lines are fits to the first order 
exponential decay. 



